A new method for producing ultracold molecular ions 
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We propose a new method for the production of ultracold molecular ions. This method utilizes 
sympathetic cooling due to the strong collisions between appropriately chosen molecular ions and 
laser-cooled neutral atoms to realize ultracold, internal ground-state molecular ions. In contrast to 
other experiments producing cold molecular ions, our proposed method efficiently cools both the 
internal and external molecular ion degrees of freedom. The availability of an ultracold, absolute 
ground-state sample of molecular ions would have broad impact to fields as diverse as quantum 
chemistry, astrophysics, and fundamental physics; and may lead to the development of a robust, 
scalable quantum computer. 



I. INTRODUCTION 



To date, the complex internal structure of molecules 
has precluded the application of laser cooling, the work- 
horse of atomic physics, as a means to produce cold 
molecular samples. Thus, several techniques have been 
developed to provide cold molecules. The most promis- 
ing of these techniques are Stark deceleration buffer 
gas cooling , and photo-association of ultracold atoms 
[3|. Despite enormous efforts on these fronts, success 
has been limited. The direct cooling methods - Stark 
deceleration and buffer gas cooling ~ are very general 
techniques that produce cold molecular samples from a 
large class of molecules, but their usefulness is hampered 
by the relatively warm temperature (^100 mK) and low 
density (~10^ cm~'^) they produce. Likewise, because 
photo-associated molecules are necessarily comprised of 
atoms that can be laser-cooled, photo-association is not 
amenable to many of the interesting precision measure- 
ment and quantum chemistry experiments. 

An exciting alternative to these methods is the study 
of cold molecular ions. Because molecular ions are easily 
trapped for many minutes, there are possibilities for cool- 
ing and interrogation that are not applicable to neutral 
molecules. Furthermore, because ion trapping is essen- 
tially species independent, many of the major goals of 
cold polar molecule physics can be realized with ions. 
The most noteworthy of these goals include: the study 
of cold chemistry, which not only has important impli- 
cations for understanding the formation of interstellar 
clouds 3, but will allow for investigation, and possi- 
ble control, of reactive collisions in the truly quantum 
regime 0); precision measurement of molecular transi- 
tions, which can be used to very sensitively measure par- 
ity violating effects [1], as well as constrain the possi- 
ble variation of the fundamental constants Q; and, per- 
haps most interestingly, the implementation of a scalable, 
quantum computation architecture 0]. Molecular ions 
possess many of the advantages that molecules have for 
quantum computation with the added benefit of simple, 
reliable trapping. A recent proposal Q has shown how 
neutral polar molecules trapped near the surface of a su- 
perconducting microwave strip-line cavity would realize 



an excellent, scalable quantum computation architecture. 
This proposed quantum computation scheme relies on the 
coupling of the molecular rotation to the electromagnetic 
field of a superconducting strip-line resonator. If strong 
coupling can be achieved between the molecule and the 
field, then it was demonstrated that both the implemen- 
tation of the necessary gates for quantum computation 
and the construction of quantum memory are feasi- 
ble for this scalable molccule-mesoscopic system. The 
attractive quantum computation scheme of Ref. Q may 
also be applicable to cold molecular ions [l^, which can 
be easily produced at the requisite density, as long as the 
strong-coupling regime can be reached. 

To date, cooling of molecular ions has been demon- 
strated through sympathetic cooling collisions with both 
cold helium buffer gas [ll| and laser-cooled atomic ions 
Unfortunately, neither of these methods can simul- 
taneously produce both ultracold and internal ground- 
state molecules - i.e. the ground state in both the in- 
ternal and external degrees of freedom, a prerequisite for 
the aforementioned experiments. While loading of molec- 
ular ions into traditional ion traps at cryogenic temper- 
atures has been demonstrated with a He buffer gas, this 
technique, which produces limited cooling of the inter- 
nal molecular degrees of freedom [l^, e.g. rotation 
and vibration, results in relatively warm molecules (T ^ 
4 K). Conversely, sympathetic cooling of molecular ions 
via co-trapped, laser-cooled atomic ions produces no in- 
ternal state relaxation due to the long-range nature of the 
Coulomb interaction [Tsj . but quickly cools the internal 
degrees of freedom down to the ultracold temperatures 
possible with laser cooling. While it may be possible 
to follow He buffer gas loading with sympathetic cooling 
via laser-cooled atomic ions to produce ultracold, internal 
ground state molecular ions, the method is unattractive 
due to the poor efficiency of He collisions for relaxing 
the internal degrees of freedom. (In Ref. [H, [13 it is 
shown that only 1 collision in 10'* to 10® collisions leads 
to any internal state relaxation.) We are aware of one 
other well-developed proposal [3, [l3| that addresses the 
need to cool the molecular internal degrees of freedom; 
however, this method is comparatively very complicated 
to implement, as it relies on a combination of Raman 
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lasers, spontaneous emission, and black-body radiation 
to producing rotational cooling, and has yet to be demon- 
strated. 

Here we propose a new method, based on sympathetic 
cooling via collisions with ultracold neutral atoms that 
should allow simultaneous cooling of both the external 
and internal molecular ion degrees of freedom. To our 
knowledge, this method of cooling has not been previ- 
ously proposed. This is likely due to the fact that the 
ionization potential, IP, of atoms that are amenable to 
laser cooling is far less than the energy released by do- 
nating that electron to the majority of positively charged 
molecular ions. Therefore, when most molecular ions col- 
lide with a laser-cooled neutral atom, charge exchange 
occurs and the molecular ion may become a sig- 
nificantly more energetic, neutral molecule. Nonethe- 
less, as detailed in this work, many cases exist for col- 
lisions between neutral atoms and certain molecular ions 
where charge exchange is energetically forbidden and the 
neutral-ion collision can only lead to cooling of the exter- 
nal and internal molecular ion degrees of freedom. Ref. 
[isl l has proposed investigating neutral atoms collisions 
with both atomic and molecular ions; however, in the 
cases considered charge exchange was possible and the 
relaxation of molecular internal degrees of freedom was 
not discussed. 

In what follows, we detail the possibilities for sympa- 
thetic cooling of molecular ions via laser-cooled neutral 
atoms. We provide calculations of the expected cooling 
rate for a specific molecular ion-neutral atom combina- 
tion, and give a brief description of a simple experimental 
design for investigating these neutral-ion collisions. 



II. METHODS 

The crux of our proposed cooling method is identifying 
combinations of molecular ions and neutral laser-cooled 
atoms that, due to energy conservation, cannot undergo 
charge exchange during a collision. Towards this goal. 
Tab. U shows a partial list of the IP of several neutral 
atoms and molecules. Since the energy released when a 
positively charged molecular ion accepts an electron is 
equal to the IP of the formed neutral molecule, charge 
exchange is energetically forbidden if the neutral atom IP 
is larger than the neutral molecule IP. Note that the num- 
bers reported in Tab. U are in eV, therefore the kinetic 
energy of even room temperature atoms and molecules is 
negligible on this scale. 

As shown in Tab. |T1 there are many possibilities for 
sympathetic cooling of molecular ions by laser-cooled 
neutral atoms. In principle, every molecular ion consid- 
ered here could be sympathetically cooled by ultracold Be 
atoms. However, laser cooling of neutral Be has not yet 
been demonstrated, due to the inconvenient wavelengths 
and low room temperature vapor pressure. For the sake 
of this proposal we choose to work with Yb atoms, which 
represent a compromise between high ionization poten- 



tial and experimental practicality, ultracold Yb atoms 
are routinely produced via laser-cooling with standard 
diode laser technology and since they are easily ionized 
to produce Yb+ , they offer an in situ method of ion trap 
diagnosis. All ions preceding Ti2 in Tab. I are energeti- 
cally forbidden to undergo charge exchange with neutral 
Yb and are thus amenable to our cooling method. 

Any two-body chemical reactions between these ions 
and neutral Yb, other than the already precluded charge 
exchange, would necessarily involve either removing an 
additional electron from the molecular ion or be an 
exchange-type reaction. In the former case, to donate 
an electron from the molecular ion would require the 
liberation of more energy than the IP of the molecular 
ion. Since double ionization potentials are typically large 
(> 10 eV) and the electron affinity of Yb is only ^^0.5 eV, 
this is forbidden for every combination known to us. In 
the latter case, reactions of the form AX+ + B BX+ 
4- A may occur if the dissociation energy of the result- 
ing BX+ ion is larger than the AX+ ion. Molecular ion 
dissociation energies are not well measured, and the pos- 
sibility of this reaction must be evaluated on a case by 
case basis; however, since the B atom has been chosen for 
its large ionization potential and this reaction necessarily 
involves removing an electron from the B atom, this reac- 
tion is usually precluded. Furthermore, by choosing the 
cation and ultracold atom to be the same, i.e. AX"'" -I- A 
AX+ -I- A, the effects of this mechanism can be mit- 
igated. Finally, energetically allowed addition-type reac- 
tions may exist, e.g. Ca -I- CaF"'" Ca2F"'" is exothermic 
by ~2.81 eV [22] , however, they require a 3-body coUision 
to conserve both energy and momentum and thus are not 
a concern at the densities considered in this work. 

From our survey of available molecules, we choose to 
work with BaCl"*" in this proposal for the following rea- 
sons: 

• BaCl"'" should be easily produced via ablation of a 
BaCl2 target. 

• BaCl^ possesses a simple ro-vibronic structure. 

• The dissociation energy of BaCl"*" is 4.8 eV, while 
the dissociation energy of YbCl"*" is 4.1 cV, thus 
the exchange reaction is precluded [2l|. Further- 
more, the mass of BaCl"*" is equal to that of neutral 
Yb and therefore momentum exchange during the 
sympathetic cooling process is favorable. 

We expound each of these points in turn below: 

A. BaCl+ Production: 

We expect that laser ablation, as detailed in Ref. [1^ 
for the production of CaF"*" from CaF2, will efficiently 
produce BaCl"*" from BaCl2. (In fact, CaCl"'" ions pro- 
duced via laser ablation have recently been trapped [2J|.) 
BaCl2 sputtering targets are readily available from many 
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Atom 


IP [eV] 


A [nm] 


Molecule IP 


[eV] 


Molecule IP 


[eV] 


Molecule IP 


[eV] 


Molecule 


IP [eV] 


Molecule 


IP [eV] 


Be 


9.3 


235 


Rba 


3.5 


Li2 


5.1 


YH2 


6.2 


TiS 


7.1 


Nal 


7.64 


Mg 


7.6 


285 


CS2 


3.6 


SrF 


5.3 


HoF 


6.2 


GeCl 


7.2 


RbBr 


7.7 


Yb 


6.3 


399 


Ka 


4.1 


SrCl 


5.5 


Ti2 


6.3 


Csl 


7.25 


Pd2 


7.7 


Ca 


6.1 


423 


KNa 


4.4 


SrBr 


5.5 


us 


6.3 


SiF 


7.26 


C6H7N 


7.7 


Sr 


5.7 


461 


KLi 


4.6 


Sri 


5.5 


ErF 


6.3 


GeBr 


7.3 


(Aniline) 




Li 


5.4 


671 


Bal 


4.7 


CaBr 


5.6 


TiO 


6.4 


SnF 


7.4 


LiD 


7.7 


Na 


5.1 


591 


BaBr 


4.8 


UO 


5.7 


Va 


6.4 


SnBr 


7.4 


CsBr 


7.72 


K 


4.3 


771 


SrF 


4.9 


CaF 


5.8 


CaO 


6.5 


Mn2 


7.4 


PbBr 


7.8 








CeO 


4.9 


CaH 


5.9 


SnCl 


6.6 


Si2 


7.4 


Bsi 


7.8 








BaF 


4.9 


CaCl 


6 


TiO 


6.8 


MgCl 


7.5 


MgF 


7.8 








PrO 


4.9 


DyF 


6 


BaO 


6.9 


PbF 


7.5 


N(CH3)3 


7.82 








Na2 


4.9 


TaO 


6 


UN 


7 


PbCl 


7.5 


(Trimethylamine 










LaU 


4.9 


im 


6 


mb 


7 


JN (01120113)3 


7.5 




7.85 








BaCl 


5.0 


SrO 


6.1 


ZrH2 


7 


(Triethylamine) 




CU2 


7.89 








LiNa 


5.0 


ZrO 


6.1 


InSe 


7.1 


HfO 


7.5 


PtSi 


7.9 








VO 


5.0 


Cal 


6.1 


TiS 


7.1 


GeF 


7.5 


Ge2 


7.9 








NdO 


5.0 


UC 


6.2 


Rbl 


7.1 


InTe 


7.6 


ZrN 


7.9 



TABLE I: Partial list of the ionization potential of candidate neutral atoms and molecules. The first three columns are data 
for atoms amenable to laser cooling. The third column, labeled A [nm], lists the required wavelength for laser cooling. The 
remaining columns list candidate molecules and their ionization potential, IP. Since the molecular IP represents the amount 
of energy released when the positively charged ion of that same molecule accepts an electron, charge exchange is forbidden in 
collisions with a neutral atom if the molecule IP is smaller than the IP of the neutral atom. We emphasize that this is only a 
partial list of molecules and that the number of molecular ions that can be cooled by this technique is likely much larger. Data 
taken from [H, H^, HJ] . 



suppliers. InRef. [2^ ablation was performed on an opti- 
cal quality sample of CaF2 for the investigation of optical 
damage via ultra-violet laser light. We envision follow- 
ing the work of Ref . [2^ where ablation of a crushed and 
pressed sample was used to improve ablation efficiency. 
One advantage of this technique is that Yb (or other Yb 
precursors, e.g. YbCl2) can be added to the sample be- 
fore it is pressed. Thus, Yb+ ions will also be produced 
in the ablation process. This will make it possible to 
load Yb"*" ions into the ion trap, allowing for trap di- 
agnosis, molecular ion-atomic ion sympathetic cooling, 
and the possibility of an efficient detection scheme of 
trapped molecular ions through fluorescence detection of 
the atomic ions [HH^,!!!. 



B. BaCl+ Structure: 

The positively-charged alkaline earth monohalides ions 
(MX+) have similar electronic structure to alkaline earth 
chalcogens, e.g. BeO, and are an ionically-bonded com- 
plex of two closed-shell atoms (M^+X~) - this property 
is both the reason that they exhibit reduced reactivity 
compared to other molecular ions and that their neutral 
atom precursor have such low IPs. Their ground-state 
is "'^E+ and thus they are described by Hund's coupling 
case (a). By analogy to the alkaline earth chalcogens 
the first excited state is though depending on the 

details of each molecule there may exist a lower-lying 
■^11+ state. Though the alkaline earth monohalides have 
been studied more than most molecular ions, very lit- 
tle spectroscopic data is available for any member. The 



case of BaCl"*" is somewhat aided by the observation from 
pyrotechnics that BaCl+ can be used to produce green 
colors in fireworks. From this observation it is likely that 
the BaCl+ A^J2^^^ ^X^S^yj transition wavelength lies 
within the range between 505 nm and 533 nm [29l |. As 
will be presented in a later section, knowledge of only the 
BaCl^ ro-vibrational structure, not the electronic struc- 
ture, is necessary to our experiments. Therefore, accu- 
rate spectroscopic identification is not necessary. The 
BaCl+ ro-vibrational spectra can be estimated from the 
known constants for BaCl, where the vibrational split- 
ting is w ^ 8.4 THz [2^ and the rotational constant is B 
^ 1.25 GHz js^. Since ionization of an alkaline earth 
monohalide results in a more tightly bound molecule, 
these values topically increase by ~10% (from a compar- 
ison of Ref. [3l| and [l^l) for the molecular ion over the 
neutral species. Nonetheless, we use these values in the 
remainder of this work as order of magnitude estimates. 

C. Sympathetic Cooling: 

1. External degrees of freedom 

The temperature, T, of a molecular ion undergoing 
sympathetic cooling with an ultracold Yb buffer gas at 
temperature TBuffer, depends on the number of colli- 
sions, n, as [3^ : 

T = TBuffer + [Tlnitial ~ Tbu/ fer)e~"'^ (1) 

with Tlnitial the initial molecular ion temperature and 
K = ^'^tiT''^ , where M is the mass of the molecular ion 
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Number of Collisions 

FIG. 1: (color online) Number of collisions required to ther- 
malize the ion motion with laser-cooled Yb (m = 173 amu) 
atoms. The shaded regions shows the range of temperature 
vs. number of collisions for all of the molecular ions con- 
sidered here that can be cooled via laser-cooled Yb atoms; 
from the best-matched BaCl'*' (m = 173 amu) to the worst- 
matched Li^ (m — 14 amu). The green shaded region rep- 
resents cooling with Yb at 1 mK, produced via cooling on 
the ^Pi ^^So line. The gray shaded region represents cool- 
ing with Yb atoms at T= 10 piK produced from second stage 
cooling on the intercombination line ^Pi <— ^Sq. 



and m is the mass of the Yb atom. A plot of Eq. 1 is 
shown in Fig. [1] for the range of possible temperature 
versus collision number for all of the molecules consid- 
ered in Tab. [T] that are amenable to sympathetic cooling 
via ultracold Yb; from the best-matched BaCl+ (m = 
173 amu) to the worst-matched LiJ (m = 14 amu). We 
have displayed two curves in this graph, one (green) for 
collisions with Yb atoms at 1 mK, produced via cooling 
on the ^Pi ^""^So line, and the other (gray) for Yb atoms 
produced from second stage cooling on the "^Pi ^^Sq 
intercombination line at 10 fxK. 

In this plot, we have conservatively used T initial = 
1000 K for the initial temperature of the ions, which is 
the approximate temperature of species created via laser 
ablation; however, as detailed in a later section, we en- 
vision first loading the molecular ion trap with a cryo- 
genic Ne buffer gas before sympathetically cooling the 
ions with the ultracold atoms. In this scheme, the ini- 
tial ion temperature would be only ^^30 K, reducing the 
number of collisions required for thermalization by a fac- 
tor of ^2. Though the cooling is optimal for BaCl"'", as 
will be shown explicitly later in this section, the number 
of required collision for all considered molecular ions is 
experimentally feasible. 

We now estimate the total elastic scattering rate as a 
function of temperature. The Eikonal approximation is 
valid when the interaction potential does not vary con- 
siderably over the scale of the de Broglie wavelength of 
the scattering particle. This approximation is appropri- 
ate as long as the collision energy, is much greater 
than the interaction strength - that is, as long as the 
idea of a semi-classical path is reasonable and the wave- 



function can be replaced by a semi-classical wavefunc- 
tion. With the application of the optical theorem, the 
Eikonal approximation gives the total elastic scattering 
cross-section as [ssj : 

<JE^k = ^Im{f{k,k)), (2) 

with the imaginary part of the forward-scattered ampli- 
tude given by: 

/•OO 

lTn{f{k,k)) ^ k 6(1 - cos(2A(6)))d&, (3) 
Jo 

where the phase of the semi-classical wave is found in 
terms of the impact parameter, 6, as: 



m --^j ^ v{V^T^)dz. (4) 

In general, a molecular ion produces a polarization in a 
nearby neutral atom, leading to an interaction potential: 

where a is the atomic polarizability, e is the absolute 
value of the electron's charge, and r is the atom-ion sep- 
aration. With straight-forward integration of Eqs. (3) 
and (4), we find the total elastic cross-section, as pre- 
dicted by the Eikonal approximation, to be 

where T is the gamma function [IJ] and v can be approx- 
imated by the molecular ion velocity. Finally, we have 
the total elastic scattering rate constant, K, as: 

where () denotes thermal averaging and fi is the reduced 
mass. The weak T^/^ dependence of the scattering rate 
constant means that sympathetic cooling is efficient even 
to the lowest temperatures. 

A log-log plot of the scattering rate constant is shown 
in Fig. [2] for the case of BaCl+ colliding with Yb, where 
/i = 173/2 amu and a = (47reo)179 a^ [3a|, for the range of 
temperatures of concern to our project. From this figure 
we see that the scattering rate constant is very large by 
atomic physics standards (a result of the relative strength 
of the ion monopole field) and for typical ultracold atomic 
densities oi p ^ 10^^ cm^"^ one would expect total scat- 
tering rates of 7 = 100 Hz - 100 kHz. Thus, the external 
degrees of freedom of the molecular ion are conservatively 
expected to sympathetically cool to ~100 jiK within only 
a few hundred milliseconds. Since molecular ions can be 
trapped for hours, this would be a very effective cooling 
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FIG. 2: Log-log plot of the elastic scattering rate constant for 
BaCl^ + Yb collisions as given by the Eikonal approximation. 



mechanism. Furthermore, the vibrational state energies 
of typical ion traps are on the order of 10 MHz [s^, thus 
the temperatures represented by this cooling mechanism 
correspond to the molecular ions primarily occupying the 
absolute ground state of the external trapping potential. 

An estimate of the range of validity of the Eikonal ap- 
proximation can be made by defining the length scale of 
the interaction potential, rchan as equal to the value of r 
where the ion-neutral interaction is equal to the incoming 
kinetic energy: 



' Char 



(8) 



The requirement that the dc Broglie wavelength be less 
than this characteristic length scale, leads to a constraint 
on the temperature for the validity of the Eikonal approx- 
imation as: 



T > 



SOOAiK. 



(9) 



Thus, the estimates made in this section should be accu- 
rate over the majority of the range of interest. 



2. Internal degrees of freedom 

As aforementioned, neither collisions with a helium 
buffer gas nor collisions with laser-cooled atomic ions 
are efficient at cooling the internal degrees of freedom 
of a molecular ion. In the former case, the inefficiency is 
due to the small electric polarizability of the He atom, 
which prevents it from coupling the molecular ion ro- 
vibrational states in the latter case, the Coulomb 

repulsion between the molecular and atomic ions results 
in long-range collisions that do not couple the molecu- 
lar ion ro- vibrational states • Neutral atomic Yb can 
scatter at close range from a molecular ion and possesses 
a polarizability ^200 times larger than that of atomic 
He. Since the interaction potential is linear in the polar- 
izability, from Fermi's Golden Rule we expect that Yb 



may be ~ 10^ times more efficient than He at quench- 
ing the molecular ion internal degrees of freedom. Given 
that data for collisions of molecular ions with He shows 
that only 1 in 10"* to 10^ collisions relax the molecular in- 
ternal states, while more polarizable buffer gases lead to 
higher quenching rates we expect that almost ev- 

ery short-ranged collision with Yb should result in some 
relaxation of the molecular internal state. 

We now calculate the expected short-range (inelastic) 
collision rate between BaCl+ and Yb. For two colliding 
particles, the energy-dependent scattering cross-section 
for the £th partial wave, with projection m, from state i 
to state / in all outgoing partial waves f , m! is generally 
written as: 

'Ji,m{E,i^ f)^^Y. \Ti,n.,i',m'iE,t^ f)\^. (10) 
i',m' 

Here Te^m,e',m'iE,i —>■ f) is the 'T-matrix', whose el- 
ements represent the probability amplitude for a tran- 
sition from the incoming spherical wave, ^'i.£.rn, to the 

out-going wave, ^'/,£',rn', and k = \J is the magni- 
tude of the wave-vector at collision energy E. Since we 
are concerned with inelastic collisions we must sum over 
all final states /, in addition to the normal sum over £, m. 
Furthermore, if we assume that any collision that pene- 
trates to short-range is inelastic with unit probability, we 
can rewrite Eq. [10] as [s^l : 

a{E, i)^Yl * ^ /) = E ^(2/ + 1)Pt[E, £), 

f.l.m I 

(11) 

where Pt^E^I) is simply the probability of transmission 
to short-range. 

Since the molecular ion-neutral atom interaction is 
spherically symmetric (neglecting the much weaker 
dipole terms) we have the 1-D Schrdingcr equation as: 



h'^i{i + i) C4 



-^]^ = E^. (12) 



2fi dr^ \ 2/ir2 

Rescaling the problem with x = r/(3, where /3 = 
\/2^Ci/h , and defining = /{2^0^) we have: 

(f^ ft{l+l) 1 \ E ^ 

- — + ( - ) * = — (13) 



E4 



In this dimensionless form Eq. (4) represents the 'univer- 
sal' Schrdinger equation for a particle with dimensionless 
energy E/E4 interacting with an potential. Pt{E, £) 
is calculated by numerically solving the Schrdingcr equa- 
tion for the potential and assuming that any flux that 
does not reflect from the potential is transmitted to short- 
range and then completely lost to inelastic processes. 
The 'universal' solution of this problem, in terms of the 
scaling factors, is given in Fig. 3. To ensure convergence 
over the calculated range we have included up to £ = 8. 
The individual partial waves are shown in upper panel 
of this graph and may be identified by their values at 
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FIG. 3: Universal scattering rate constant [K/ 0^) vs. Energy 
{E/E4). Upper panels shows individual partial wave contri- 
butions, labeled by its respective £ value, for up to ^ = 8. The 
lower panel shows the total scattering rate constant, with the 
inset depicting the observance of the Wigner threshold scal- 
ing law for low temperature inelastic collisions. The scatter- 
ing rate for a specific pair of colliding ions and atoms can be 
found by calculating the /3 parameter, which depends only on 
the reduced mass and atomic polarizability. 



E/E4 = 2000, where the scattering rate constant is high- 
est for £ ^ 8 and decreases monotonically to the lowest 
value for £ = 0. Consistent with the Wigner threshold 
scaling law for low temperature inelastic scattering, our 
calculations find a oc 1/k. As a result, the scattering rate 
constant becomes constant as the energy goes to zero, as 
seen in the inset on the lower panel of Fig. [31 

As shown in Fig. 3, the inelastic scattering rate is 
roughly constant over the entire energy range considered 
and has a value of ~ 6 x 10^ /P^. For BaCl+ collisions 
with Yb, where fi = 173/2 amu and a = (47reo)179 a^ 
[3^ . this yields a temperature independent scatter- 
ing rate constant of A' = 4 x 10~^ cm'^s"^. For a typical 
ultracold atomic density of p = 10^"'^ cm~'^, one would 
then expect a total inelastic scattering rate of ^100 kHz. 
Thus, in only a matter of milliseconds it should be pos- 
sible to quench any internal excitation of the molecules 
into their absolute ground state. Furthermore, even if 
the assumption of unit efficiency internal relaxation for 
Yb collisions is off by orders of magnitude the scheme is 
still viable ~ i.e. if only 1 in 100 collisions produced re- 
laxation K Ri 1 kHz. With these rates it may be possible 



to scatter multiple photons from a single molecular ion 
since the inelastic collisions will relax the molecular ion 
ro-vibrational states to the ground state, serving as an 
effective dark-state repumper. Interestingly, this could 
lead to both efficient fluorescence detection and, possi- 
bly, direct laser cooling of the molecular ions. 

As a check of the model, we predict a quenching rate 
constant for collisions between He and BaCl+ of K ~ 
10"^'^ to 10""'^^ cm'^s^^. This is consistent with values 
reported in the literature for similar collisions [l3| . 

Finally, to cool a sample of diatomic molecular ions 
from T « 1000 K to ~1 mK, requires removing an 
amount of energy E = 5/2NMkBT, where Nm is the 
total number of molecular ions and fcs is the Boltzmann 
constant. Approximating the ultracold Yb atoms as an 
ideal gas we expect the temperature of the Yb atoms to 
rise by an amount AT = 5/3Nm/NaT, where A'^ is the 
total number of trapped Yb atoms. For typical values of 
Nm « 10^ and Na « 10^ AT « 1 mK. Thus, the heat 
capacity of a typical ultracold atomic sample is enough to 
cool the molecular ions down to ultracold temperatures. 
Furthermore, if the estimates here are not accurate, e.g. 
the ions are produced at a higher temperature, the cool- 
ing process, as described in the next section, could simply 
be repeated with a fresh 'batch' of ultracold Yb atoms. 



III. EXPERIMENTAL IMPLEMENTATION 

Though charge exchange is energetically forbidden for 
collisions between ground state (^Sq) Yb and BaCl+, Yb 
atoms in the excited ^Pi state are not forbidden from 
donating an electron to the molecular ion. Therefore, 
all sympathetic cooling collisions between the neutral 
atoms and molecular ions must occur with the Yb atoms 
in their ground state. This can be arranged by using 
a 'dual' version of the linear ion trap described in Ref. 
[3^. As shown in Fig. IDJa), this trap design consists of 
four linear quadrupole electrodes in the standard config- 
uration; the radio-frequency trapping voltage is applied 
to the two diametrically-opposed bare electrodes, while 
the other two electrodes are grounded. Three pairs of 
short DC electrodes are slipped over the two grounded 
quadrupole rods. The DC electrodes are separated from 
the grounded rods by a thin insulator. By controlling the 
DC voltage applied to these 'end-cap' electrodes, ions can 
be trapped in either the left or right 'chamber', as well 
as shuffled between the two. 

In a manner similar to that described in Ref. (39j 
molecular ions, produced via laser ablation, can be loaded 
into the left chamber of the ion trap (represented by 
the red cloud in Fig. [H^a)). The ion loading efficiency 
may be improved over that of Ref. (39| by using a cryo- 
genic Ne buffer gas [i^ at T « 30 K. Once the ion trap 
is loaded, the cryostat temperature may be lowered to 
'^A K, resulting in fast cryopumping of the Ne buffer gas 
by charcoal sorbs [4l[. Next, ultracold Yb atoms are 
collected in a standard MOT (represented by the blue 
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cloud in Fig. [UJa)) and then transferred into an opti- 
cal dipole trap whose axis is aligned with the ion trap 
axis (the blue cloud in Fig. IHb)). Once the ultracold 
atoms are collected, the molecular ions are transferred 
from the 'collection' chamber into the 'collision' chamber 
by control of the voltage applied to the DC electrodes. By 
simply grounding the center DC electrode the ions will 
move to the right in Fig. [3) Once the ions have moved 
to the rightmost chamber, reapplying the DC voltage to 
the center electrode results in a ion trap overlapped with 
the ultracold atomic Yb (represented by the overlapping 
blue and red clouds in Fig. HJc)). Achieving neutral 
atom densities of 10^'^ cm~^ over the entire ion trap vol- 
ume is routine with an optical dipole trap [i^l and thus 
the estimates of cooling rates from the previous section 
are conservative. 

The result of a trajectory simulation of the cooling pro- 
cess for a single BaCl+ ion is shown in the bottom panel 
of Fig. 01 Here the ion 'temperature' [i^ is plotted ver- 
sus time. In the first shaded region, corresponding to Fig. 
mja), the molecular ion is sympathetically cooled from 
the laser ablation temperature (conservatively assumed 
to be 10^ K) down to the Ne buffer gas temperature. 
Approximately, 100 fis later the buffer gas is removed 
from the simulation and the transfer process begins (Fig. 
m^b)). Initially the ion gains kinetic energy, but by timing 
the transfer process appropriately, it can be reconverted 
into potential energy. Once the transfer process is com- 
plete, the molecular ions come into equilibrium with the 
10 fiK Yb buffer gas. These trajectory simulations in- 
clude the ion micro-motion and model the sympathetic 
cooling collisions as randomly directed hard-sphere colli- 
sions. For collisions with the Ne buffer gas a conservative 
cross-section of 2x 10~^^ cm^ is used, while the Yb cross- 
section is calculated from Eq. [6] 

To carefully study the BaCl'^/Yb coUisional properties 
it is necessary to monitor the molecular ion external and 
internal temperatures. The motional temperature of the 
molecular ions may be measured throu gh o bservation of 
the fluorescence from co-trapped Yb+ [H, [l^, [H • Af- 
ter the sympathetic cooling stage, Yb"*" ions can be cre- 
ated from ionization of the ultracold Yb atoms and the 
residual Yb atoms removed by extinguishing the opti- 
cal dipole trap. Because of the long-range nature of the 
Coulomb interaction, the Yb"*" and BaCl+ temperatures 
will quickly come to equilibrium. Thus, by simply mea- 
suring the atomic lincshape the molecular temperature 
may be inferred. This in situ thermometry will allow 
continual monitoringof the BaCl+ temperature during 
the cooling process 14J] . 

The molecular ions' internal temperature, i.e. rota- 
tional and vibrational population distribution, can be 
measured using the technique of resonantly-enhanced 
multi-photon ionization as proposed in Ref. |15| . It may 
also be possible to measure the internal state population 
with a method similar to that developed for the spec- 
troscopy of ^"^Al fisj . 




FIG. 4: (color online) (a) BaCl^ ions are trapped in the left 
chamber, while Yb atoms are captured in a MOT centered on 
the right trap chamber, (b) Yb atoms are transferred into an 
optical dipole trap and the DC electrode voltages are switched 
to transfer the ions into the right chamber, (c) The DC elec- 
trode voltages are reconfigured to trap the molecular ions in 
the rightmost chamber, where they are overlapped with the 
ultracold Yb atoms. The bottom panel shows the temper- 
ature of a BaCl^ ion during the three stages of the cooling 
process. Each shaded region, with is label (i), corresponds to 
the similarly labeled panel above. 
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IV. OUTLOOK AND SUMMARY 

Even when energetically allowed, charge exchange 
rates are typically much slower than elastic collisions 
rates [l^- Therefore, the estimates of this proposal arc 
likely very conservative and efficient sympathetic cooling 
of molecular ions may be possible with more easily pro- 
duced ultracold Rb or Cs atoms. Similarly, it may not be 
necessary to quench the ultracold atoms to their ground 
state, allowing continual cooling of the molecular ions by 
the MOT. This is especially attractive since without con- 
tinued laser cooling the RF-trap induced micro-motion 
of the molecular ions will eventually lead to sympathetic 
heating of the atoms. (Though advantageous, continued 
cooling is not necessary. Due to the large number of 
atoms, this heating is negligible on the time-scales neces- 
sary to quench the ion internal motion and could easily 
be overcome by sympathetic cooling with a co-trapped 
atomic ion.) 

Interestingly, our proposed technique naturally ex- 
tends to a combination of an 'on-chip' MOT, similar to 
that proposed in Ref. [46| . and a surface ion trap, similar 
to those proposed in Ref. Thus, using lithographic 

techniques it should be possible to construct simple and 
scalable 'on-chip' ultracold molecular ion traps that will 
provide molecular ions in their absolute internal and ex- 
ternal ground states, i.e. the requisite starting point for 
the exciting experiments discussed in the introduction. 



In summary, we have presented a new and promis- 
ing mechanism for the production of ultracold molecular 
ions. The key realization of this work is that the com- 
monly held belief that molecular ions cannot undergo pri- 
marily elastic collisions with atoms that are amenable to 
laser cooling, is not an absolute truth. Wc have shown 
that for carefully chosen ion-atom pairs, charge exchange 
is energetically forbidden and thus sympathetic cooling 
with ultracold neutral atoms is indeed possible. There- 
fore, this proposal represents a dramatically simplified 
scheme for producing ultracold molecular ions in their 
absolute internal and external ground state as compared 
to previous proposals [1^1 , which rely on a four- fold com- 
bination of sympathetic cooling with atomic ions, laser 
induced Raman transitions, spontaneous emission, and 
black-body radiation. Because ion trapping is species 
independent, many of the major goals of cold polar 
molecule physics can be realized with ions. Thus, the 
impact of this proposal is clear. With an efficient means 
of producing ultracold molecular ions, interests as di- 
verse as quantum chemistry, astrophysics, fundamental 
physics, and quantum computation may be pursued. 
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